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By tempered Monte Carlo simulations, an Almeida-Thouless (AT) phase-boundary 
line in site-diluted Ising spin systems is searched for. Spins interact only through 
dipolar fields and occupy a small fraction of lattice sites. The spin- glass susceptibil- 
ity of these systems and of the Sherrington-Kirkpatrick model are compared. The 
correlation length as a function of system size and temperature is also studied. The 
results obtained are contrary to the existence of an AT line. 
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I. INTRODUCTION 

The collective behavior of some spin systems is controlled by dipole-dipole interac- 
tions. It is so in some magnetic nanoparticle^ arrays,^ in some crystals of organometal- 
lic molecules,^ as well as in some magnetic salts, such as LiHoF 4 . In LiHoF 4 , uniaxial 
crystal-field anisotropy forces the Ho ion spins to point up or down along the anisotropy 
axisP^ A model of Ising spins with dipole-dipole interactions ought therefore to capture 
the main features of the magnetic behavior of LiHoF4. This system orders ferromagneti- 
cally at low temperatures, which, as Luttinger and Tisza 7 showed long ago, is accidental. 
Had the Ho ions crystallized in a simple cubic lattice, for instance, it would have ordered 
antiferromagnetically. 8 This illustrates how delicate the balance between dipolar fields com- 
ing from different sources is. The frustration that underlies such a balance is expected 
to lead to spin-glass behavior in disordered-Ising-dipole (DID) models which mimic the 
LiHo x Yi_ x F4 family of materials^ if x <C 1. 

Some details about LiHo :r Yi_ x F4, such as the symmetry of its crystalline lattice, are 
irrelevant if x <C 1. Other details, such as transverse fields, which have no place in 
the DID model, do make a difference. Thus, interesting quantum effects that have been 
observe d 11 * 12 ^ in LiHo x Y 1 _ x F4 at low temperatures are beyond DID models. On the other 
hand, a clear picture of the DID model seems like a good starting point for the study of 
quantum dipolar systems. Thus far, no such clear picture exists. 

Several experiment j 11 ! 13 ! on LiHo a; Y 1 _ a .F 4 suggest there is a paramagnetic (PM) to spin 
glass (SG) phase transition when x < 0.25, but some skepticism remains. 12 Some computer 
simultion of DID models 1 ^ point to a PM phase for all nonvanishing temperatures. However, 
the opposite conclusion has been drawn more recently) 10 * 15 ! 

Below the transition, the nature of the hypothetical SG phase of DID models remains 
rather unexplored. Simulations for zero applied field suggest!^ the DID model behaves in 
three dimensions (3D) somewhat similarly to the XY model in 2D. Thus, 3 would be the 
value of the lower critical dimension di of DID models in zero applied field. Note, however, 
(i) that the correlation length of the Edwards-Anderson (EA) model has previously been 
observed^ to behave similarly, as a function of system size and temperature, (ii) that cLl < 3 
was nevertheless drawn from this behavior, and that (iii) this fits in with a 4 - 2.5 value 
that has recently been inferred for the EA model from other evidence.^^ I know of no 
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reported work on the behavior of DID models under applied longitudinal magnetic fields. 

Whether there is a thermal phase transition, between the PM and SG phases, as the 
temperature T is lowered in an applied magnetic field H is an important question. An 
H — T phase-boundary line was long ago discovered in the Sherrington-Kirkpatrick 21 (SK) 
model by de Almeida and Thouless (AT).^l For its existence in the EA model, there is 
both favorable^ 2 ^ and unfavorable^ 2 ^ evidence. To get a feeling for the physics involved, 
consider first the argument of Fisher and Husej^S which in turn follows from Imry and Ma's 
argument)^ for the instability of diluted AFs to the application of a magnetic field. Consider 
a patch of n spins in a SG state at H = 0. Because all the nearest neighbor bonds are of 
random sign, the numbers of spins pointing in opposite directions are then expected to differ 
by ~ n 1 / 2 . The Zeeman energy therefore changes by AEh ~ Hn 1 / 2 if a patch of n spins is 
flipped when H 7^ 0. Let the corresponding energy change coming from broken bonds be 
given by AEj ~ Tn / 3 , which defines the stiffness^ constant T and the stiffness exponent 0. 
Fisher and Huse^ further showed that < (d—l)/2 for the EA model (more recent numerical 
work gives^l 1.21n(0.4d) for 1 < d < 6), whence AEj < AE H follows for a sufficiently 
large value of n. Widespread spin reversals of this sort on macroscopic systems would lead 
to a state with a q = overlap with the initial state. (The spin overlap q between two spin 
configurations may be defined as the total fraction of sites on which spins point in the same 
direction minus the fraction of sites on which spins point oppositely.) Because dipole-dipole 
interactions are long ranged, the above argument is not immediately applicable to the DID 
model. Data for the mean square q2 of the overlap between equilibrium states at H = 
and atpS H = 0.2 is exhibited in Fig. [I^l for the DID model, for x = 0.35, all T and various 
system sizes in 3D. These results suggest that indeed q — > as L — > 00 for the DID model 
as well. Analogous results are shown in Fig. [T]d for the SK model. Again, q — > as L — > 00 
seems to ensue. This is in spite of the fact that an AT line is known to exist for the SK 
model. Whereas Imry and MaP^ could conclude that a small magnetic field can destroy the 
antiferromagnetic phase of a dilute antiferromagnet (AF), the analogous conclusion could 
only be drawn for the DID model if it were known to fit the droplet scenaricP^] (j n which 
there is no ground state degeneracy). This is why Fig. [T^i provides insufficient evidence for 
the nonexistence of the Almeida-Thouless line in the DID model. An analogy with a simpler 
system is helpful at this point. 

Consider an isotropic AF. Upon the application of an arbitrarily small magnetic field H, 
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FIG. 1. (Color online) (a) Plots of c/2 vs T/x for an L 3 DID model on a SC lattice, for x = 0.35 
and the L values shown in the graph. Here, q2 is for two replicas, both of which are in equilibrium 
but under different applied fields, Hi = and H2 = 0.2. (b) Plots of q2 vs T/J for the L 3 SK 
model, the L values shown in the graph, and Hi = and H2 = 0.4 for replicas 1 and 2. In both 
(a) and (b), most error bars do not show because they hide behind icons. Lines are guides to the 
eye. 

all spins rotate uniformly till they point nearly perpendicularly to H. From a canted AF 
alignment, spins can better minimize the ground state energy. It takes a nonvanishing H to 
further drive this "spin-flop" phase beyond the H-T boundary line, into the paramagnetic 
phased This is illustrated in Fig. |2]for the XY model in 3D. The phase transition on the H- 
T boundary line can take place because the applied field does not completely lift the ground 
state degeneracy. Two degenerate states (for two sublattices) survive. Fluctuations between 
these two states enable the existence of an H-T boundary line. Thus, sublattice symmetry 
is broken below the H-T line. Analogously, critical fluctuations between various low energy 
states take place on an AT line. In the SG phase, different replicas of a SK system can 
stay on different states. This sort of replica equality breaking, is known as replica symmetry 
breakin(j^ (though no symmetry operation relates these states). 

In the droplet scenario there are only two states, related by global spin inversion. An 
arbitrarily small magnetic field therefore lifts this degeneracy. Only one state survives, 



5 




FIG. 2. (Color online) H/J vs T/J, where J is a nearest neighbor exchange constant, phase 
diagram of the antiferromagnetic XY model in 3D. Data points come from MC simulations. The 
fitting curve, H = 11.8(1 - T/T c ) a37 , where T c = 2.20J, is also shown. 



which leaves no room for critical fluctuations to occur at any nonzero H. Hence, Fisher and 
Hus d 29 * 34 -* concluded that q = between two states, one at H = and another one at H ^ 0, 
implies the state for H ^ is not a SG state. Thus, the nonexistence of an AT line is a 
clear cut prediction of the droplet model. 

The aim of this paper is to establish whether there is an AT phase-boundary line in a 
site diluted DID model in 3D. This is to be done by means of the tempered Monte Carlo 
(MC) method.^ The behavior of the DID model, has been previously showrP^ to depend on 
x and T only through T/x for s<l. It therefore suffices to study how the model behaves 
as a function of T and if at a single value of x. 



A brief outline of the paper follows. The DID model is defined in Sec. II A The boundary 



conditions are described in Sec. II B The definition of the spin-overlap parameter q and how 



it is calculated can also be found in Sec. II B How equilibration times of the DID model 



under tempered MC rules are arrived at is described in Sec. II C Equilibrium results for 
the spin-glass susceptibilities Xsg of the DID model and SK models, both for H = and and 



H 7^ 0, are compared in Sec. Ill Equilibrium results for the correlation length £l of the 



DID model are also given in Sec. Ill Results for both Xsg and £l are clearly in accord with 
the absence of an AT phase-boundary line in the DID model. Further concluding remarks 



appear in Sec. IV 
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II. MODEL, METHOD AND EQUILIBRATION 

A. Model 

The DID model on a simple cubic (SC) lattice is next denned. All dipoles point along 
the z axis of the lattice. Each site is occupied with probability x. The Hamiltonian is given 

by, 

ij i 

where the sums are over all occupied sites, except for % = j in the double sum. Oi = ±1 on 
all occupied sites i, 

T ij = e a (a/r ij ) 3 (l-34/r? j ), (2) 

Tij is the distance between i and j sites, % is the z component of r^, e a is an energy, and 
a is the SC lattice constant. 

For H = 0, the DID model has been shownM to have an equilibrium PM-SG transition 
if x < 0.65(5) (in SC lattices). Furthermore, the PM-SG transition temperature is given by 
T sg = l.0(l)x for all x < 0.5. 

For comparison, a few results for the SK model are also shown. Then, all exchange con- 
stants are given random values chosen independently from the same Gaussian distribution, 
centered on with a J 2 /N mean square deviation. 

Unless otherwise stated, all temperatures and energies for the DID model are given in 
terms of e a /kB and e a , respectively. The magnetic field H is defined by Eq. ([TJ to be an 
energy, and is therefore also given in terms of e a . All times are given in MC sweeps (MCS). 

B. Method 

I use periodic boundary conditions (PBC), in which a periodic arrangement of replicas 
span all space beyond the system of interest. These replicas are exact copies of the Hamil- 
tonian and of the spin configuration of the system of interest. A spin on site i interacts 
through dipolar fields with all spins within an L x L x L cube centered on it. No interac- 
tions with spins beyond this cube are taken into account. (Additional details of the PBC 
scheme used here can be found in Ref. El) This may seem odd, because dipolar interactions 
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make themselves felt over macroscopic distances. That is why different "demagnetization 
factors" apply to differently shaped macroscopic bodies.^ On the other hand, demagnetiza- 
tion factors vary with system shape, but not with macroscopic system size. Indeed, the error 
that is introduced by this method was shown in Ref. [TDlto vanish as L — » oo, regardless of 
whether the system is in the paramagnetic, AF or SG phase (but not near a ferromagnetic 
phase transition). This enables us to disregard interactions of any one spin on site i with 
any spin beyond an L x L x L cubic box centered on site i. 

In order to bypass energy barriers that can trap a system's state at low temperatures 
the parallel tempered MC algorithm is used herej^ following the steps outlined in Ref. [TUl 
Configuration swap rates between systems at temperatures T and T + AT were checked to 
be reasonably large throughout. The smallest swap rates ensue for the lowest temperature 
(i.e., T = 0.05) and the largest systems (i.e., L = 10). Then, swap rates in equilibrium were 
found to be approximately 0.3, i.e., 30% of all attempts made for configuration exchanges 
are successful. Swap rates increase slowly with increasing T in the spin-glass phase, and 
faster above T sg . 

In order to be able to calculate spin overlaps between different equilibrium states at the 
same temperature, not one, but two sets, each one of n identical systems, are allowed to 
evolve independently in parallel. All In systems start from independently chosen random 
configurations. The temperature spacing AT between systems in each set was chosen to 
be AT = 0.05. Checks for equilibrium are described below, following the time dependent 
spin-overlap definitions. 

As usual, the Edwards- Anderson overlap^ between identical systems (replicas) 1 and 2 
is defined by, 

q = N- l ^v (3) 



where 



■ S = °?°f\ (4) 



crj and cr| are the spins on site j of identical replicas (1) and (2) of the system of interest. 
Unless otherwise stated, identical replicas have, as usual, the same Hamiltonian. Exception- 
ally, for Figs. [T^i and [T]d, different fields Hi and H 2 are assumed to be applied to replicas 1 
and 2, respectively. 
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C. Equilibration 

The purpose of this subsection is to establish how long it takes the DID model to come to 
thermal equilibrium. In order to be able to follow the equilibration process (under tempered 
MC rules), some useful quantities are next defined. First, two replicas are allowed to evolve 
independently, starting at t = from two uncorrelated random states and x v . Let 
Q2{t | t^jty) be the average of q 2 at time t over all sample realizations. Different samples 
start from different random pairs of states, t M and x v . In q 2 (t \ t^,t u ), t M and x v appear only 
to remind us that all initial pairs of states at t — are uncorrelated random states. 

During equilibration, q 2 (t \ t^,t u ) is expected to increase up to its equilibrium value, q 2 . 
In Fig. gfi, q 2 (t | t M ,t„) is given for T/x = 0.571 and T/x = 1.14, at H = 0. In Fig. Qy, 
H = 0.2, but everything else is as in Fig. [3^. 

Finally, assume two replicas start evolving independently from the same equilibrium state 
c M at time t = 0. That is, any state c M is selected from the sequence of states the system 
of interest goes through after thermal equilibrium has been reached. The time dependent 
equilibrium correlation function q 2 (t | e^, e^) is the average of q 2 at time t over all sam- 
ple realizations. Again, e M , e M appear in q 2 {t | e M , e M ) only to remind us that both replica 
evolutions start at t — from the same e M equilibrium state. 

Note that q 2 (0 \ e M , e M ) = 1, and that ergodicity implies q 2 (t \ e^, e^) — > q 2 as t — > oo. 
Therefore, q 2 {t \ e M , e M ) is expected to be an upper bound to q 2 . Plots of q 2 (t \ t^,c^) are 
shown in Fig. ||i for T/x = 0.571 and T/x = 1.14 at H = 0. In Fig. go, H = 0.2, but 
everything else is as in Fig. [3^,. 

A measure r q of equilibration times in tempered MC evolutions, under the conditions 
specified in Table I, is defined graphically in Fig. [3^i. It turns out that r q 10 2 , 3 x 
10 3 , 5 x 10 4 , 10 6 for L — 4, 6, 8, 10, respectively, for the DID model. For equilibrium 
observations below, all MC runs went on for 2r s MCS. Values of t s are given in Table I. They 
fulfill r s ^> r q . Equilibrium was achieved in the first half of each run, that is while t < r s . 
All time averages for the calculation of equilibrium values were taken while r s < t < 2t s . 

The following rules for the time evolution of q 2 (t | r M , x v ) under a tempered MC algorithm 
are noted in passing. The first rule, q 2 (Q \ v^t,,) = 1/N, which follows from the fact that 
spin configurations are initially random, is exact. The second rule, that q 2 (t \ t^t u ) ~ fa( T ' 
when 10 < t < r q , and (l(T) ~ 0.4 (weakly dependent on T and L), follows from plots 
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FIG. 3. (Color online) (a) Plots of q2(t n \ Vfj,,x v ) and of q 2 (t | e M ,e^) vs t, in MC sweeps, for the 
values of T/x shown, x = 0.35, L = 8, and H = 0.0. The procedure that was followed to arrive at 
values for r q is illustrated. For qiit n | e^, e^), equilibration was allowed to proceed for 5 x 10 5 MC 
sweeps before measurements were taken. Lines are guides to the eye. Error bars are given by the 
size of the icons, (b) Same as in (a) but for H = 0.2. 

of q 2 {0 | t/i) tv) vs t, such as the ones shown in Figs. [3^, and [3Jd. Further digression into 
equilibration behavior under tempered MC rules is beyond our aim here, which is simply to 
determine equilibration times. 



III. EQUILIBRIUM RESULTS 

Equilibrium results obtained from tempered Monte Carlo simulations are reported in this 
section. These results are for both site-diluted DID models and SK models. The SK model, 
in which an AT line is known to exist, is examined for comparison purposes. 

All the data given here for DID models is for x = 0.35. This is well below x c (~ 0.65), 
in a regime where DID models on SC lattices have been showrP^ to have an SG phase if 
H = 0. Furthermore,^ T sg = 1.0(l)x. 

In the search for the existence of an AT line in DID models, I apply well known criteria.^ 



10 



TABLE I. Number t s of Monte Carlo sweeps (MCS) allowed, first for equilibration and, subse- 
quently, for averaging over equilibrium, and number N s of samples for the SK model and for DID 
models of various linear sizes L. For the SK model, 1? is the number of spins. For the DID model, 
L is given in units of the lattice constant, each site is occupied with 0.35 probability, the tempera- 
ture T fulfills 0.05 <T< 2.0, and the temperature spacing between systems in the tempered MC 
runs is AT = 0.05. 



Model 


L 




N s 


SK 


4 


500 


10 3 


SK 


6 


10 3 


10 3 


SK 


8 


5 x 10 3 


10 3 


SK 


12 


10 4 


500 


DID 


4 


10 4 


5 x 10 3 


DID 


6 


10 5 


3 x 10 3 


DID 


8 


10 6 


10 3 


DID 


10 


10 7 


300 




Let 



ij 



(5) 
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FIG. 5. (Color online) (a) Plots of Xsg/L 2 vs T for the DID model, for x = 0.35, H = 0, and the 
L values shown in the graph, (b) Same as in (a) but for H = 0.1. Error bars show wherever they 
protrude beyond icons. 

where 50j = 0, — (0,), and k = (2ir/L, 0, 0), perpendicular to all spin directions. Note x(0) 
is the spin-glass susceptibility, Xsg- 

In the paramagnetic phase, short range spin-spin correlations imply Xsg is finite, but 
Xsg — > oo as the PM-SG critical point is approached. At the critical point, Xsg/L remains 
finite as L — > oo in the SK modelP^ Plots of Xsg/ L vs T, shown in Fig. [4^, for H = and 
various values of L, exhibit the expected behavior. Similar plots for H = 0.4 are shown 
in Fig. [3|d. Clearly, Xsg/L curves for various values of L do cross, as expected for the SK 
model, at a nonvanishing value of T. Furthermore, they do so at T/J = 0.55(5), which is, 
within errors, on the AT line J 22 ' 44 1 

For DID models, one must first decide how to scale Xsg- Recall that, quite generally, finite 
size scaling predicts a finite limit of Xsg/L 2 ~ v at a critical point as L — > oo. Furthermore,^ 
7] ~ in DID models. Plots of Xsg/L 2 vs T for H = and various values of L are seen to 
cross, as expected, at T/x ~ 1.0 in Fig. |5^i. 

Not knowing in advance the value of rj for the hypothetical AT line in DID systems, 
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universality is next assumed. Thus 77 = is assumed to hold for H 7^ as well. To probe 
for an AT line, I vary T with H > constant. One does not want to miss the AT line by 
choosing too large a value of H. I let H = 0.1. Since T sg ~ x for x < 0.5 and H = 0, 
and x = 0.35 has been chosen everywhere, H — 0.1 gives a Zeeman energy of Q.3ksT sg 
approximately, which is a rather small field. (For comparison, recall that H along the AT 
line increases beyoncP^ H = 3ksT sg as T — > in the SK model.) 

Plots of Xsg/L 2 vs T at H = 0.1 are shown in Fig. [5|d. These results show the AT line, 
if there is one, is restricted to H < 0.1, that is, to | H |< 0.3T sg . 

If instead of r/ = one uses rj = —0.3, from the table given in Ref. [3B]for the EA model in 
3D, the plots in Figs. |5^i and|5jD are slightly modified. For H = 0, curves for different values 
of L would then cross at T/x = 0.8, instead of at T/x = 1.0, as in Fig. [5ji. For H = 0.1 the 
main effect is to spread all curves shown in Fig. [5]d further apart, thus strengthening the 
conclusion drawn above about the AT line. 

The correlation length £ is more convenient than Xsg to work with, because £/L — > 
remains finite at the critical point as L — > 00 while £/L — > in the paramagnetic phase. 
Diagnostics with £/L is thus free from errors in the value of rj. Let 

where k is a unit vector along k, and the L subscript is a reminder of the fact that, in- 
evitably, the sum in the equation is performed over finite L 3 size systems. Obviously, £ L is 
a correlation length measured along the k direction. 

Numerical computations of the double sum in Eq. ^ are however time consuming. In 
addition, £ L is not well defined if (5<f)i5<j)j) decays (as it doeskin the SG phase) more slowly 
than r~ p and p < d + 2. Both difficulties are avoided with the definition,^ 



a- 1 



4sin 2 (fc/2) 



x(Q) 

X(k) 



(7) 



Note that £l — > £ L / \/2 as £l/L — > in the macroscopic limit if £l is finite, since (i) exp ik • 
can then be replaced by 1 + ik ■ ry — (k ■ rjj) 2 /2 in Eq. ([5]), and (ii) (2/fc) sin(/c/2) — > 1 then. 
Thus, Eqs. ^ and ([7]) are qualitatively equal in the paramagnetic phase. Equation ^ is 
therefore, as has become customary in SG workp^^ adopted here as the definition of 
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o I , , , , , , , , i I 

1 T/x 2 

FIG. 6. (Color online) (a) Plots of i L /L vs T for the DID model, for x = 0.35, H = 0, and the L 
values shown in the graph. Error bars show wherever they protrude beyond icons, (b) Same as in 
(a) but for H = 0.1. (c) Same as in (b) but for ^ instead of (,l/L. 

correlation length. 

In the paramagnetic phase, C,l/L — > as L — > oo. What various assumptions about the 
SG phase imply for the variation of ^l/L with L is discussed in some detail in Sec. VB of 
Ref. HQ] In short, (i) < 3 (recall d^ is the lower critical dimension) implies C,l/L — > oo 
(and a nonvanishing Xsg/N) in the SG phase as L — > oo, (ii) d^ = 3 implies ^l/L remains 
finite (and Xsg/N — > but Xsg oo) in the SG phase as L — > oo. 

Plots of £,l/L vs T for the DID model at H = and x = 0.35 are shown in Fig. [6^. The 
L — > oo limit of £l/L, obtained from 1/L — >■ extrapolations of £l/L in Ref. QUI is also 
shown in Fig. [6^i. 

Plots of £,l/L vs T for the DID model at H = 0.1, are shown in Fig. ^jp. Curves do 
cross for the smaller values of L, but the trend is reversed for larger L. Then, ^l/L decrease 
as L increases, at least for the temperatures studied. With a confidence level above 99%, 
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95%, and 85%, £ L /L is smaller for L = 10 than for L = 8, at T/x > 0.43, T/x = 0.28, and, 
T/x = 0.14, respectively. As for x sg above, this is the behavior one expects of if there 
is no AT line. 

Plots of £x vs T on Fig. [6J3 are perhaps more revealing. Clearly, £^ saturates for all T to 
a finite value for L > 8, as one expects from a paramagnetic phase. 

We end this section with a comment about spatial anisotropy in DID systems. Recall 
interactions along the z-direction, parallel to the spins axis, are twice as large as in a per- 
pendicular direction. The "longitudinal" (for k along the ^-direction) correlation length £ l L 
is consequently somewhat larger, up to twice as large for high temperatures, than the trans- 
verse correlation length £l. More importantly, £l/L appears to suffer from finite size scaling 
corrections in a way that £l/L does not: whereas £l/L for systems of various sizes all cross 
at approximately the same temperature in Fig. [6^1, £ l L /L do not quite do so for L = 4, 6, 8 
and 10. The crossing points for £ l L /L drift towards T sg as system sizes increase. For this 
reason, transverse correlation lengths are more convenient to work with. For H > 0.1, we 
find no qualitative difference between £} L and £l. 

IV. CONCLUSIONS 

Spin-glass behavior in an applied magnetic field H has been studied. More specifically, 
I have numerically probed a site-diluted Ising dipole model of magnetic dipoles for the 
existence of an Almeida-Thouless phase-boundary line. This DID model has been previously 
shown 10 to have, in three dimensions, at H = and low temperatures, (i) an AF phase for 
x > x c , where x c = 0.65(5), (ii) a (marginal) SG phase for x < x c , that is cLl — 3, (iii) a 
behavior for 1 Ci c that is independent of lattice structure and depends on x and T only 
through T/x, and (iv) T sg /x ~ 1. Spin-glass behavior as a function of T and H can therefore 
be inferred for all x <^ x c from that at a single small value of x. 

Here, equilibrium results, from tempered Monte Carlo simulations, are reported for 
Xsg/L 2 ~ v and ^l/L for the DID model at x — 0.35, various temperatures and system sizes, 
at H = and H = 0.1. The criterion for the existence of an AT line, that Xsg/L 2 ~ v and 
£l/L be independent of L at the critical point, has been shown here to work well for (i) the 
SK model at H = and H = 0.4J, that is, H = 0AT sg , for which the answer has long been 
knownP as well as (ii) for the DID model at H = 0. For H = 0.1, that is, H ~ 0.3T S9 , the 
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trend observed in the data is clearly away from Xsg/L 2 ^ v or £l/L becoming independent of 
L at any T as L — > oo. Indeed, £ L saturates to a finite value beyond L ~ 8 for all T. All of 
this is consistent with the absence of an AT phase boundary line in the DID model, at least 
above H > 0.3T C . 
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